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Abstract Mixing processes such as convection, over- 
shooting and rotational mixing have long been known 
to affect the evolutionary properties of low-mass stars. 
While modeling a 1.2 M star, we encountered a semi- 
convective region outside the fully convective core, rem- 
iniscent of the well-known situation for massive stars. 
In this study we focus on low-mass stars presenting con- 
vective cores and, by applying different prescriptions 
for the determination of the convective boundaries and 
using different mixing descriptions for the dynamical 
processes in the core, we look for the effects of semicon- 
vective mixing in the interior structure of the stars and 
its observable quantities. With this purpose, we have 
constructed different sets of evolutionary models using 
a stellar evolution code (GARSTEC), and analyzed the 
models looking for imprints of these processes. 
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1 Introduction 

It has long been known that different mixing processes 
in stellar cores strongly influence the overall evolution of 
stars. The treatment of convectively unstable zones and 
overshooting regions in evolutionary models affects the 
outcoming luminosity, temperature, and main-sequence 
lifetime of the stars, among other consequences. 
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One of the mixing processes that has been extensively 
studied in massive stars evolution is that of semicon- 
vection. During the hydrogen-burning phase the ra- 
diation pressure makes the core expand. Opacities 
are increased outside the core by electron scattering, 
where a chemical discontinuity appears as a result 
of this core expansion. Sin ce the pioneer work of 



Schwarzschild fc Harm! ([19581 ) , several authors have in- 
vestigated the occurrence of semiconvective mixing in 
massive stars and its effects on stellar evolution (e.g. 
Stotherdfl975 IStothers fe Chm1ll975t lLangerlll985h . 
Although semiconvection was initially thought to oc- 
cur only in massive stars, it was also found in low-mass 
stars as a consequence of a discontinuity in the molec- 
ular weight produced either by convective core expan- 
sion due to the increa sing importan c e of the C NO cy- 
cle over the pp chain (jMitalas! Il972c I Said Il975l ) , or by 
the retreating co nvective core leaving behi nd a chemi- 
cal discontinuity ( Faulkner fc Cannonll973lh both cases 
producing higher opacities outside the convective core. 
We address the issue of the determination of convective 
boundaries and the treatment of mixing zones which 
present a gradient in the molecular weight by the inclu- 
sion of the Ledoux criterion for convective instability, 
and a diffusive approach for semiconvective mixing in 
our stellar evolution code. 



2 Mixing Prescriptions 

The convective zone boundaries are defined either with 
the Schwarzschild criterion (Schwarzschild fc Harm 



1958 ) or the Ledoux criterion ( Ledoux 19471) . For the 
latter case, the Ledoux temperature gradient is defined 

as 
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where V p =tZln^/rflnP, and (3 is the ratio of gas 
pressure to total pressure. If the Ledoux criterion is 
applied, a zone is considered convectively unstable if 
Vl < V ra( j, while in the case of Schwarzschild crite- 
rion this reduces to V a( j < V ra( j. Convective zones can 
be mixed instantaneously or as a diffusive process us- 
ing the convective velocity estimated from the mixing- 
length theory. 

The inclusion of molecular weight gradients in the cri- 
terion for convection allows us to identify semicon- 
vective zones, which are defined as the regions where 
V a( j < V ra d < Vl- Once a semiconvective region has 
been found, the temperatur e gradient is calcula ted us- 
ing the prescription given bv lLanger dl98ill985h . and a 
time-dependent diffusive mixing in the semiconvective 
zone is performed. Both the diffusion coefficient and the 
temperature gradient in this prescription depend on an 
efficiency parameter of semiconvection (a sc ). We stress 
the fact that unlike this study, previo us efforts to treat 
semic onvection in low- mass stars fe.g. lCrowe fe Mitalas 
1982) considered the mixing in the semiconvective layer 
to be performed by adjusting the composition until con- 
vective neutrality according to the Schwarzschild crite- 
rion is reached, making the hydrogen discontinuity in 
the boundary of the convective core disappear. 
Convective overshooting is also considered in our study 
and it is implemented as a diffusive process consisting in 
an exponential decay of the convective velocities within 
the radiative zone. The diffusion constant is given by 
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where £ corresponds to an efficiency parameter cal- 
ibrated with open clusters, H p is the pressure scale 
height, z is the distance from the convective border, 
and the constant Do is derived from MLT-convective 
velocities. The extent of the overshooting region is lim- 
ited in the case of small convective cores. 



(|Angulo et al.lll999t) . We computed models without 
diffusion and a He content of Y= 0.25. 
Several models were computed to test the different 
mixing prescriptions and the convective boundary def- 
inition. Within a convective zone the mixing is per- 
formed instantaneously, while in a semiconvective zone 
or overshooting region the mixing is carried out as ex- 
plained in Sect. [2J Calculations were made for both the 
Schwarzschild and the Ledoux criterion for the defini- 
tion of the convective zones, with and without including 
extra mixing due to overshooting and semiconvection. 
We explore the effects of this processes in models start- 
ing from the pre-main-sequence phase and evolved up 
to hydrogen exhaustion in the core. 
In Fig. [TJ we present the convective core evolution dur- 
ing the main-sequence lifetime for three models of 1.2, 
1.5 and 2.0 M Q . The left panel shows the case where 
we applied the Schwarzschild criterion for convective 
boundaries determination, while the right panel shows 
the case where the Ledoux criterion was applied. In 
both cases, solid lines depict the size of the convec- 
tive core, while the dotted lines represent the extent of 
the semiconvective region. No mixing was performed 
throughout the semiconvective zones (a sc — 0), which 
do appear in models for both convective criteria. We 
choose these models as they represent cases of inter- 
est: one model where the development of the convec- 
tive core depends on the convective boundary defini- 
tion used (1.2 Mq), one model with growing convective 
core during the main-sequence lifetime (1.5 Mq), and 
one model with a receding convective core during the 
H-burning phase (2.0 Mq). However, it is important 
to keep in mind that the effects of semiconvection are 
washed out when overshooting is included as described 
in Sect. [H at least for the calibrated value of the effi- 
ciency parameter £, no matter which criterion is used 
to define the boundary of the convective regions. For 
the cases where it is present, the convective core dis- 
appears when the central hydrogen is exhausted, which 
corresponds to the end of the main-sequence phase. 



3 Models 

For our model calculations we used the Garching Stell ar 
Evolution Code CGARSTEC. IWeiss fe Schlattl (ho08h . 
The input phy sics considered incl udes the OPAL equa- 
tion of state ( Rogers et al.lll996 ) complemented with 
the MHD equation of st ate for low temperatures 
( Hummer fe Mihalasl Il988l). Ferguson' s opacities for 
low temperature ([Ferguson et al.l 120051) . OPAL opac - 
itie s for high tempera t ures (llglesias fc Rogers! Il996l ). 
the iGrevesse fc Sauvall (|1998l ) solar mixture, and the 
NACRE compilation for thermonuclear reaction rates 



In this study, we will focus mainly on the 1.5 Mq case, 
as it corresponds to the model which has the most ex- 
tended semiconvective zone of all the considered ones. 
In Fig. [2] we present the main-sequence evolutionary 
models for this mass value, including now semiconvec- 
tive mixing. As can be seen, already with a very small 
value for the efficiency parameter of semiconvection, 
this extra mixing makes the convective core grow much 
closer to the size of the Schwarzschild criterion case. 
This result is expected, as large values of a sc imply 
faster mixing, and in the limit when a sc — > oo the mix- 
ing is performed instantaneously and the Schwarzschild 
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Fig. 1 Left panel: Models computed using the 
Schwarzschild criterion for convective boundary definition. 
Right panel: Models computed using the Ledoux criterion 
for convective boundary definition. Solid lines depict the 
convective core evolution, while dotted lines show the ex- 
tent of semiconvective zones. No semiconvective mixing or 
overshooting was applied. See text for details 

limit should be recovered. As was already mentioned, 
we consider the main-sequence lifetime to end when hy- 
drogen is exhausted in the center, which coincides with 
the center of the star becoming radiative. It is inter- 
esting to note that in our models, the lifetime on the 
main-sequence is reduced when the Ledoux criterion 
and semiconvective mixing are applied to the models 
with respect to the usual Schwarzschild criterion, which 
is a n opposite result with r espect to previous findings 
(e.g. Crowe fc Mitalaslll982 ). This is because in the pre- 
vious implementations of semiconvection for low-mass 
stars, the mixing through that region recovering con- 
vective neutrality gave a fresh input of hydrogen fuel 
supply to the core and extended the H-burning phase. 
This is not the case when a diffusive, time-dependent 
mixing is considered for the semiconvective region as in 
our models. 

In Fig. [3]we present a Hertzsprung- Russell diagram for 
the 1.5 Mq models considered. Marked with diamonds 
are models which have very similar positions in this di- 
agram, but the mixing prescriptions applied make their 
evolutionary stage in the main-sequence different. This 
is shown in Fig. |4j where the hydrogen profiles are 
plotted for the selected models. The sharp chemical 
composition gradients are clearly seen, signature of the 
location of the convective core boundary. It is inter- 
esting to note that for the model computed with the 
Ledoux criterion and no mixing in the semiconvection 
region, the location in the HRD corresponds to a star 



Fig. 2 Convective core evolution for models of 1.5 Mq. 
Black lines present the cases of Ledoux criterion with no 
semiconvective mixing (dots) and Schwarzschild criterion 
(solid). Blue lines show the models where additional mix- 
ing was included: semiconvection with a sc — 0.001 (dots), 
a BC = 1 (dashed), and overshooting (solid). 



almost at the end of its main-sequence lifetime, which 
is not the case for the other models considered being in 
a less advanced evolutionary stage. 



4 Conclusions and Future Prospects 

In the present study, we have explored the effects of dif- 
ferent convective boundary determinations and mixing 
prescriptions on the internal structure and evolution of 
low-mass stars. We have done so considering mainly 
the case of a 1.5 Mq star, for which we have shown the 
difference in the internal structure resulting as a conse- 
quence of the applied physics, although the position in 
the Hertzsprung-Russel diagram remains very sim ilar. 

we 



In an upcoming study (|Silva Aguirre et al.l [201oh 



will further investigate the effects of these processes for 
different masses and metallicities, as well as changes in 
the efficiency parameters of semiconvective mixing and 
overshooting, and the inclusion of diffusive mixing. It 
is also important to keep in mind that semiconvection 
will play an important role in the evolution of a star de- 
pending on the transition between the pp chain and the 
CNO cycle as the main source of energy production (for 
a given metallicity) . This transition is critically depen- 
dent on the value of the (N 14 + p) cross section, which 
has been s ubject of substant ial reductions in the past 
years (e.g. iMarta et al.ll2008t ). 
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Fig. 3 Hertzspmng-Russell diagram for the 1.5 M© mod- 
els. Black lines present the cases of Ledoux criterion with no 
semiconvective mixing (dots), and Schwarzschild criterion 
(solid). Blue lines show the models where additional mix- 
ing was included, semiconvection with a BC = 0.001 (dots), 
a sc = 1 (dashed), and overshooting (solid). 



Although accurate stellar parameters are of course of 
vital importance, asteroseismology can also be used as 
a tool to disentangle the observational degeneracy by fo- 
cusing on the structural differences among the models 
produced by the mixing proc esses applied. It has been 
show n that both p-mod es dPopielski fc Dziembowski 
20051 ) and g-modes (e.g. iMiglio et al.l l2008h are sensi- 
tive to changes of the molecular weight in the interior 
of stars, and could in principle be applied to distinguish 
the type of mixing that has occurred in the stellar in- 
terior. 
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Fig. 4 Hydrogen profiles for the 1.5 Mq models with 
Xc~0.25, marked with diamonds in Fig. [3] Black lines 
present the cases of Ledoux criterion with no semiconvec- 
tive mixing (dots), and Schwarzschild criterion (solid). Blue 
lines show the models where additional mixing was included, 
semiconvection with a sc = 0.001 (dots), a sc = 1 (dashed), 
and overshooting (solid). 
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